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Abstract 

Physical decomposition of the non-Abehan gauge field has recently solved the two-decade-lasting 
problem of a meaningful gluon spin. Here we extend this approach to gravity and attack the 
century-lasting problem of a meaningful gravitational energy. The metric is unambiguously sep- 
arated into a pure geometric term which contributes null curvature tensor, and a physical term 
which represents the true gravitational effect and always vanishes in a flat space-time. By this 
decomposition the conventional pseudo-tensors of the gravitational stress-energy are easily rescued 
to produce definite physical result. Our decomposition applies to any symmetric tensor, and has 
interesting relation to the transverse-traceless (TT) decomposition discussed by Arnowitt, Deser 
and Misner, and by York. 

PACS numbers: 11.15.-q, 04.20.Cv 



1 



Gauge invariance is the most elegant and efficient principle for constructing interactions 
in the present field theories of physics. By requiring field equations to be gauge invariant, 
the manner of the couplings (and self-couplings) of various fields are strongly constrained. 
This applies both to the standard model of the strong and electro-weak interactions, and 
to Einstein's gravitational theory. For the latter case gauge invariance refers to general 
covariance under arbitrary coordinate transformation. It is rather annoying, however, that 
a theory built uniquely out of the gauge- invariance requirement does not seem to guarantee 
gauge invariance for all physical quantities. In hadron physics, e.g., in the two-decade 
efforts to understand how the nucleon spin originates from the spin and orbital motion of its 
quark and gluon constituents, one encounters severe difficulty in finding a gauge-invariant 
description of the gluon spin and orbital angular momentum. Only recently, a solution 
was obtained in Ref. and further developed in Ref. |2|. A more celebrated and still 
unsolved gauge-dependence problem is the energy density of the gravitational field. After 
countless attempts of nearly a century, a convincing solution is still lacking. A refiection of 
this desperation is the often heard argument that, since the effect of gravity at any point 
can be eliminated by transiting to a free-fall frame, gravitational energy is intrinsically 
non-localizable and can at best be quasi-local to a closed two-surface js], Q| . 

The key obstacle to constructing all physical quantities gauge-invariantly is the inevitable 
involvement of the gauge or gravitational field together with their ordinary derivatives, which 
are all intrinsically gauge dependent. The idea in Refs. is to decompose the gauge 

field: Afj^ = Afj^ + Afj_. The aim is that will be a physical term which is gauge- covariant 
and always vanishes in the vacuum, and ^4^ is a pure-gauge term which solely carries the 
gauge freedom and has no essential physical effects (particularly, it does not contribute to 
the electric or magnetic field strength). Equipped with the separate A^ and A^, a naively 
gauge-dependent quantity (such as the gluon spin S = E x A) can easily be rescued to be 
gauge-covariant, simply by replacing A^ with A^, and by replacing the ordinary derivative 
with the pure-gauge covariant derivative constructed with instead of A^. 

Mathematically, a well-defined separation A^ = A^^ + A^^ means an unambiguous pre- 
scription for constructing A^ and A^ out of a given A^. The properties (especially, gauge 
transformations) of and A^ are then inherently determined via their mathematical ex- 



pressions in terms of An. In Refs. 



2], it was found that and can indeed be solved 



in terms of An by setting up proper differential equations and boundary conditions, which 
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lead to unique solutions for and with desired physical properties. In this paper, we 
show that this method can be generalized to gravitational theory. The metric tensor g^^^y 
is unambiguously decomposed into the sum of a physical term g^j^i,, which represents the 
true gravitational effect, and a pure geometric term g^^, which represents the spurious grav- 
itational effect associated with coordinate choice. Gauge-dependence of the gravitational 
energy originates exactly from the fact that the metric may contain a spurious gravitational 
effect. While in a flat space-time the Cartesian coordinate with vanishing affine connection 
seems a natural choice, in an intrinsically curved space-time no coordinate is obviously more 
natural than others, hence it is no longer a trivial task to get rid of the spurious gravita- 
tional effect. In an accompanying paper 5[], we discuss a gauge-fixing approach, by defining 
a unique physical coordinate which contains no spurious gravitational effect. In this pa- 
per, we present the more general field-decomposition approach, by seeking a prescription to 
identify the geometric ^^j^ for a given metric g^^^ in any coordinate. 

As for gauge theories, we find that the prescription is again a set of defining differential 
equations, which are displayed most concisely in the form [6|: 

= d,K. - d.K, + K^:. - KJ:, = 0, (la) 
9'% = 0. (lb) 

The notations require some caution: fj^^ is the purely geometric part of the affine connection. 
Its relation to (jfj^u is analogous to that of T^^ and gfj^u- 

= Ir^'id^g.u + d,g^, - d^-g,,). (2) 

Here g^^ is defined as the inverse of g^^, i.e., g^'^gpu = The aim of this choice is that R^'^^i, 
in f lTa|) is just the Riemann curvature of g^y. It must then be noted that g^^ = g^'^ — g'^" is 
not the inverse of g^^. (In fact, the physical term g^^, may not have an inverse at all.) The 
difference — f = r^^^, is defined as the physical connection. It is not related to gfj,^ as 
inEq. 

To comprehend how Eq. ([1]) is chosen, how it gives solution for T^^ and T^^ [and further 
for gfj_y and g^,^ by Eq. ([2])] with desired properties, and how the solution in turn is employed 
to solve the gauge- dependence problem of the gravitational energy, it is most helpful to 



recall the parallel constructions for gauge theories in Refs. 



2|. In Abelian case, the gauge 



field transforms as — > A' = — d^^u, which leaves the field strength invariant: 



F^u = d^Ay — dyA^ — > F'^^ = F^y. The defining equations for the separation A^ = A^ + A^ 
are 

F^u = d^A, - d,A^ = 0, (3a) 
diAi = 0. (3b) 

Eq. fl3a|) has very clear physical meaning: the pure-gauge term A^ gives null field strength. 
Eq. (]3bp can be regarded as the transverse condition for a physical photon with zero mass. 
But to avoid confusion with the radiation gauge condition diAi = for the full Ai, it is more 
helpful to think in a mathematical way that Eq. ([3]) are the needed differential equations 
to solve A^ and A^. Since + = A^, it suffices to examine A^. To this end we rewrite 
Eq. (I3al) as 

d^A, - d,A^ = F^,. dSH) 

A clever way to solve is to act on both sides with di, set fi = i and sum over i, and use Eq. 
(l3b|) . This gives 

V^i, = diF,,, or A, = -i-diF,,, (4) 

where we have required a natural boundary condition that, for a finite system, the physical 
term vanish at infinity, as does the field strength F^j,. 3] The explicit solution in Eq. 
dl]) indicates clearly that the physical field is gauge invariant, and hence the pure-gauge 
field A^ = A^ — Afj_ carries all the gauge freedom and transforms in the same manner as 
does the full A^^. Moreover, Eq. tells us that the physical term vanishes if the field 
strength F^^, = 0. 

In non-Abelian case, the gauge transformation is more complicated: A'^ = UA^U'^ — 
^UdyU'^ . The field strength now contains a self-interaction term, and transforms covariantly 
instead of invariantly: Ffj_y = dfj,Ai, — d^Afj, + ig[A^, Ay] — )■ F'^^ = UF^^W . It is fairly non- 
trivial to choose proper defining equations for the non-Abelian A^ and A^. They were 
originally proposed in Ref. and further developed in Ref. j2| to be: 

Ff,u = d^Ay - d^A^ + ig[A^, A^] = 0, (5a) 
■DiAi = diAi + tg[Ai, Ai] = 0. (5b) 

We will shortly show that Eq. ([5]) gives solution for and A^ with desired gauge- 
transformation properties: 

= uA^u\ A'^ = uA^u^ - '-ud.uK (6) 



By these properties, T>^ = (9^ + igi-A^, ] is a pure-gauge covariant derivative for the ad- 
joint representation, and Eq. (|5]) is covariant under non-Abehan gauge transformations. 
Analogous to the Abehan case, Eq. (E]) says that is a pure- gauge field giving null field 
strength, and the physical field satisfies a "covariant transverse condition". However, as 
we remarked in the Abelian case, the real justification for Eq. is that they are the right 
mathematic equations to solve and in terms of Afj_, with desired gauge transforma- 
tions in (|6]). Again, we examine with trivial boundary condition, and rewrite Eq. (|5]): 



d^A, - d,A^ = F^, + ig{[A^ - A^, A,] - [A^, A,]) (7a) 
diA, = ig[A,Ai]. (7b) 

Due to non-linearity, these are not easy to solve. To proceed, we employ the usual technique 
of perturbative expansion, which applies when either the coupling constant g or the field 
amplitude is small. For a small g, e.g., we write = Aj?-* + gA^^^ + g'^A^if^ + • • • . Eq. ([7]) 
can then be solved order by order. The zeroth-order term A^^'^ satisfy the same equations 
as fl3aip and fl3b|) . Its solution is given by Eq. (|4]), and can in turn be used to solve the 
equations for the leading non-trivial term A^f^^\ 

- = t[Af - A,, iW] - z[i(°\ A.], (8a) 

a.if) = (8b) 

The solution is obtained by the same strategy for Abelian case, and can be further employed 

" (2) 

to solve the next-order term A]i , and so on. Given validity of this perturbative expansion, 
the solution to Eq. ((Tj) is unique. This uniqueness has important implications: a) F^i, = 
is necessary and sufficient for J.^ = 0; and b) A^ and A^ have the gauge transformations as 
in (|6]). The proof of b) is as follows: Eq. ([6]) is solution of Eq. (|5]) with A^^ and A^ replaced 
by A'^ and A'^, and since the solution to Eq. (|5]) is unique, Eq. ([6]) gives the right gauge 
t r ansf ormat ions . 

We now turn to the gravitational equations ([T]) and ([2]). Because of non-linearity, we 
have to rely again on perturbative method, and require that the gravitational field be at 
most moderately strong. Namely, the magnitude of h^^, = g^^, — rj^^, (with r]^,^ the Minkowski 
metric) is smaller than 1 and can be treated as an expansion parameter. It then takes a little 
algebra to show that Eqs. ([1]) and ([2]) can be solved similar to the gauge-field equations. We 
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proceed by first looking at the physical connection f^^, to which we can assign a natural 
boundary condition that (for a finite system) f vanish at infinity as does the Riemann 
curvature -R^'g-^y We define an expansion f = f + f ^i^'' + • • • in orders of h^j,^. For the 
first-order term T^^, we get from Eq. ([1]) 

5,f^W - a.f^W = RP^^,, (9a) 
ff ) = 0. (9b) 

Solution: Set ^ = a = i in Eq. (J9al) . sum over i, and use Eq. (l9bjl . we get 

= i?^,,,. (10) 

Then, act on both sides of Eq. (|9all with c)i, set fj. = i, sum over i, and use Eq. (fTOl) . we 
obtain the solution 

f^l'^ = i(5.^''... + 5.^''..J- (11) 
This can then be employed to solve the second-order term ra^^\ From Eq. ([1]), we have 

fp(2) _ fp(2) _ /fpCi) _ -pp ^f«(l) _ fp(i)-pa _ (-tpW rp \r°'W _i_ rpWr^' (lOa) 

Here h'^^ = t]^'' — (7'^'^. Though looking tedious, Eq. (fT2|) can be solved similar to Eq. (|9]). 
The solution can be further employed to continue the perturbative procedure up to any 
desired order, in principle. 

Having separated the affine connection, we can use Eq. ([2]) to solve the metric separation, 
gf^u = (jfiu + fif^u- It is useful to define /i^;, = hf,i, + h^^, thus g^^ = Vf^u + h^,u and g^y = h^.^. 
We again look at the physical term h/j^j, which can be assigned a trivial boundary condition. 
As for f^j^, we define an expansion h^^y = h^f}} + in orders of /i^;^. From Eq. ([2]), 

we derive the first-order equation 

d;h^£ + - = 2r7,.f^« (13) 

Interchange cr, y in Eq. (fT3l) and add the result back to Eq. (fT3l) , we get 

d.U^l = ^p^VP}^^ + n,,V^}^\ (14) 

Act on both sides with di, set = i and sum over i, we obtain 

= a.(r/,.ff ) + v,.T:i'^) = Vp.R'u. + Vp.R'ua- (15) 



where in the second step we have used Eq. (fTOj) . Since this is the first-order equation, indices 
can be lowered by the Minkowski metric. Then by noticing the symmetry property of RaUv, 
we finally obtain the solution 

Here and below a comma is used to denote derivative when too many occur. The superscript 
on -R^ilj, is to remind that it is computed to first-order in h^^. Rigorously speaking, the 
second expression requires that h^^ (not just R^^^y) vanish at infinity. 
For the second-order term hau-, we derive from Eq. ([2]) 

o(2) +0^2) -0^2) =2rt VP^'^'^ + n'^fU^Mh +h -h ) 

- h^^JML + ^^alu - h^iD- (17) 

Solution of h!§}. is similar to /z.^^, though more tedious. The perturbative solution for h^y 
can be continued to the same order as F^^. 

After obtaining g^i, and g^^ = g^iy—g^y, we must remark on how gp'' and g'^" are computed. 
By definition, g^'^ is the inverse of ^^j,. Then, g^^ is computed as g^'^ — g^^- At lowest order, 
h^'^ , g'^'^, and g'^" are just related to h^^, g^^, and ^^j, by the Minkowski metric. But this 
property is lost at higher orders. 

The solutions we obtain show the desired property that the physical terms F^^ and h^i, 
vanish if and only if R''^^^, = 0, i.e., the space-time is intrinsically flat. It is also illuminating 
to look at the property of the pure geometric terms F^^ and h^^. To this end we rewrite 
Eq. ffTbj) as 5'*-'Tf^. = g^^T^j. This indicates that in order to have F^^ = (so that the 
spurious gravitational effect is absent), it is necessary that g^-'T'l- = 0. On the other hand, 
given validity of our perturbative expansion, g^^T^j = will lead uniquely to f = 0. We 
therefore name a coordinate in which g^^T^j = the "pertinent coordinate". (Similarly, in 
gauge theories, the radiation gauge diAi = leads to the solution for the pure-gauge field 
Afj^ = 0, and can be termed the "pertinent gauge" [s!].) The pertinency condition g^^T^J = 
is just what we find in Ref. js] the "true radiation gauge for gravity". It is straightforward 
to verify that the spherical coordinate is not "pertinent" even in a flat space-time. This 
explains why it gives unreasonable gravitational energy by the traditional pseudo-tensors. 
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We note that the pertinency condition is fairly non-triviaL E.g., while the Cartesian 
coordinate in flat space-time gives F^^j, = and is clearly pertinent, the quasi-Cartesian 
coordinate in a curved space-time is not necessary pertinent, e.g., the simplest Schwarzschild 
solution: ds'^ = (^ |^j|/g/2r ) ~ ~^ ^^)^ dr^. Moreover, it is not trivial to convert this 
coordinate to a pertinent one, except at linear order It is exactly the non-triviality 
of the pertinency condition that calls for our field-decomposition approach, which works 
straightforwardly in any coordinate, and can pick out the true gravitational content of the 
metric up to moderate strength. 

We are now in the position to explain how to calculate a physically meaningful energy 
density of the gravitational field, for any given g^^i, of a finite and not-too-strong gravitating 
system. The metric Qfj^i, may either be obtained by solving the Einstein equation directly, or 
may just be worked out with some guessing, or even be the experimentally measured result. 
First, the metric is put into the pertinency test: If one finds g^-^Tl- = 0, it means that this g^^^ 
contains no spurious gravitational effect, thus can be used directly in the traditional pseudo- 
tensors to compute the energy density. If, instead, g^^T^j ^ 0, it means that this g^i, does 
contain spurious gravitational effect, and one should revise a pseudo-tensor by replacing the 
quantities in it with their corresponding physical counterparts, which are obtained by the 
field-decomposition approach we just presented. This would give a concrete gravitational 
energy as physical as that in the pertinent coordinate. 

Discussion. — (i) Various pseudo-tensors show a high degeneracy concerning the total 
energy of a gravitating body. It would be interesting to examine whether such degeneracy 
persists to the level of a meaningful density. 

(ii) In gauge theories, gauge transformation and Lorentz transformation are two different 
manipulations. Therefore, in Eq. fl3bp / fl5b|) . is gauge invariant /covariant so as to make 
the equation gauge invariant /covariant. However, to make the equation hold in any Lorentz 
frame, the physical field ^4^ must not transform as a four-vector. This is an inevitable phys- 
ical feature of a massless particle with spin-1 or higher js]. In general relativity, however, 
gauge transformation and coordinate transformation mean the same thing. Therefore, to 
make Eq. fllb|) hold in any coordinate, the physical term F^^ must not transform covari- 
antly under four- dimensional transformations, even linear (Lorentz) ones. This manifests 
the masslessness of the gravitational field. But by our construction F^^ is indeed a true 
tensor under spatial transformations, following the same line as in proving the non-Abelian 
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transformations in Eq. ([6]). 

(iii) At leading order, is essentially the field defined in the "pertinent coordinate" 
as we discuss in Ref. [5|, where we have derived the second expression in Eq. f ll6p by 
a method of gauge transformation. Moreover, the expression mimics exactly the form of 
the "transverse" part of the matter stress-energy tensor, derive in Ref. [5j] by yet another 
method: 

4- = S^J - -i-{d^dkS) + d,dkS\ - d,d,S\), (18) 

where Sfj_i, = T^i, — ^ri^,^TPp. This "coincidence" is actually profound and reveals that 
our tensor-separation is a unique extension of the usual vector-separation by curl-free and 
divergence-free conditions: Riemann curvature is the unique covariant "curl" of a tensor, 
hence comes Eq. (JTall . The uniqueness of the expression in Eq. ( flbll is explained in jsl. 

(iv) Arnowitt, Deser and Misner (ADM) discussed a linear orthogonal separation of a 



symmetric spatial tensor 



10|: hij = hj^ + hfj + hfj, where hjj'^ is transverse and traceless, 



hjj is transverse, and hfj is longitudinal; all expressed uniquely via hij-. 

L 111 

= fi,j + /* = '^i.^ik,k — -^:^hki^kH) (19a) 

^5 = liS^jf - i/J), f = - (19b) 
/ij.^ = h,^-hl-hf^. (19c) 

ADM regard hfj^ as the physical part of the gravitational field. At linear order, both hfj'" 
and our hfj_u are gauge invariant. But a key difference is that in our method the rest part 
h^i, is a pure gauge, while in the ADM method hjj is also gauge invariant and only hfj is a 
pure gauge. This implies that hfj'^ does not contain all physical content of hij, and is not 
as pertinent as hfj_^. Since at linear order /i^j,, hfj'", and hjj are all gauge-invariant, we can 
expect some relations among them. Remarkably, indeed, a little algebra shows 

f = ^^ifc = hkk - i^hki^ki, (20a) 

Thus, the relation of hjj^ and hij is similar to that of the TT gauge and our pertinency 
condition: They agree for pure waves without matter source, but disagree otherwise (sf. 

(v) York has proposed a different extraction of TT component from a symmetric tensor: 
/i^.^ = hij — hij — ^Sijhkk, with hfj another longitudinal part and ^6ijhkk a trace part. [11 1 



At linear order, the explicit expression is: 

2 



K = ^i^J + ^i.* - fijWk,k (21a) 

Wi = -=r-{hik^k — -7hkk,t — -i^r-hki^kii) (21b) 
V ^ ^ V 

It can be checked that at linear order Kf^ defined by York equals that of ADM. Moreover, all 

gauge dependence is contained in the pure-gauge part Wij + Wj^i in hfj, while the — |5jjVl4,fe 

term in hfj can join ^Sijhkk to make a gauge-invariant combination: 

7T(/ifcfc — 2VFfc,fc) = -^^ihkk — ^hki^ki) = T^f^ ■ (22) 

It must be noted, however, that the Wi of York differs from the of ADM, and the 
pure-gauge terms defined by York and ADM are different: Wij + Wj^i ^ fij + fj^i. They are 
both much more complicated than our pure-gauge term in Eq. f|T6|) : 

lj,,u ~l~ ^fi ^"(^/jiji TT^ji,/^)- (23b) 



— -. {h^i^iu ~\~ hj/i^in hii^^y) (23a) 

The relations between our decomposition and that of ADM and York, especially beyond 
the linear order, will be further explored elsewhere. 12 1 
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